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ABSTRACT We review the current status of studies of large-scale structure in the X-ray Universe. 
After motivating the use X-rays for cosmological purposes, we discuss the various approaches used 
on different angular scales including X-ray background multipoles, cross-correlations of the X-ray 
background with galaxy catalogues, clustering of X-ray selected sources and small-scale fluctuations 
and anisotropies in the X-ray background. We discuss the implications of the above studies for 
the bias parameter of X-ray sources, which is likely to be moderate for X-ray selected AGN and 
the X-ray background (bx ~ 1 — 2). We finally outline how all-sky X-ray maps at hard X-rays 
and medium surveys with large sky coverage could provide important tests for the cosmological 
models. 
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1. INTRODUCTION 

In the current cosmological picture, galaxies, clusters and large-scale structures have 
grown from initial small perturbations in the density of the Universe via gravita- 
tional collapse. Cosmological models are required to meet two basic observational 
constraints: on the one hand the Universe at z ~ 1500 was very smooth, as the cos- 
mic microwave background (CMB) is seen to have anisotropies of amplitude ~ 10 -5 ; 
on the other hand local mass inhomogeneities measured through the distribution 
of galaxies exhibit fluctuations of the order ~ 1 on scales ~ 10 Mpc. Different cos- 
mologies, however, predict highly discrepant ways in which structures on different 
scales grow up to the current state from the CMB initial conditions. The largest 
discrepancies occur at redshifts z ~ 1 — 5 which is when galaxies began to collapse 
and to form stars. Accessing these intermediate redshifts will provide crucial tests 
for the cosmological models. 

The isotropy of the cosmic X-ray background (XRB) on large angular scales 
{^j- less than a few % on scales of degrees and larger) suggests that most of the X- 
ray photons we receive from the Universe must have been originated in the distant 
Universe. Surveys at different depths carried out with ROSAT have revealed that 
50-70% of the (soft) XRB is resolved into point sources, mostly Active Galactic 
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Nuclei (AGN) of different classes. Although there are still some discrepancies in the 
determination of the X-ray luminosity function and its rcdshift evolution, there is 
no doubt that most of the XRB originates at redshift z > 1. Boyle et al (1994) and 
Page et al (1996) who find their samples of X-ray selected AGN consistent with pure 
luminosity evolution models, predict a peak in the X-ray volume emissivity around 
z ~ 1.5 — 2. Miyaji et al (1998) instead find better consistency with luminosity 
dependent density evolution, in which case the X-ray volume emissivity in AGN 
more luminous than 10 44 5 ergs _1 (which for the broken power-law shape of the 
luminosity function account for most of the X-rays emitted by AGN) rises steeply 
from z = 0toz = l — 2 with no evidence for a decline at higher rcdshifts. In both 
cases it is clear that soft X-ray emission from the extragalactic sky comes mostly 
from redshifts z = 1 — 2 or larger, in a situation very similar to the star formation 
in the Universe (Madau et al 1996, Boyle & Terlevich 1998). Studying the X-ray 
Universe is then likely to provide a major handle to understand the evolution of the 
Universe at intermediate redshifts and therefore it is an issue of prime cosmological 
relevance. 

There are other reasons to prefer X-rays to carry out cosmological studies. On 
the one hand the high-latitude X-ray sky is 'clean', at least at photon energies above 
2 keV, galactic absorption has negligible effects and the contribution of the Galaxy 
to the XRB is less than a few % (Iwan et al 1982). A further reason is the small 
content in stars of high galactic latitude surveys, ranging from 25% at bright fluxes 
down to probably less than 10% at the faintest fluxes. 

In this paper we review the current status of studies of the large-scale structure 
of the Universe, which up to now has produced relevant but certainly not spectacular 
results. The two main questions that we address are: 

• Do X-ray sources (and the XRB) trace mass in the Universe and what is their 
bias parameter? 

• What are the best observational approaches to obtain information on the 
large-scale structure of the Universe at intermediate redshifts with X-rays? 

Except when otherwise stated we use Ho = 100 ftkms -1 Mpc - , qo = 0.5 and 
A = 0. 

2. THE X-RAY SKY ON THE LARGEST SCALES 

The distribution of the XRB fluctuations on the largest scales and their link to 
inhomogeneitics in the distribution of matter has been an active field of research for 
many years. The observational resources have been mostly limited to the HEAO- 
1 A2 experiment which scanned the sky with a resolution of 3° x 1.5° at photon 
energies 2-60 keV. 
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I (deg) 


b (deg) 


Err (deg) 


XRB ampl (%) 


Rcf 


CMB 


264 


48 




0.15 


F96 


AGN 


318 


38 


30 




MB90 


Clusters 


260 


5 


15 




PK98 


Soft XRB 


288 


25 


19 


1.7 


PG99 


Hard XRB 


338 


47 


25 


0.11 


S99 



TABLE 1. Dipoles of X-ray source populations and the XRB. References arc Fixsen 
et al (1996); MB80: Miyaji & Boldt (1980); PK98: Plionis & Kolokotronis (1998); 
PG99: Plionis & Georgantopoulos (1999); S99: Scharf et al (1999). 



2.1. The dipole of X-ray sources 

Since the Galaxy is moving with respect to the frame where the CMB would be 
isotropic towards I — 264°, b — 48°, there must be an overdensity of sources which 
are pulling us towards that direction. The distribution of X-ray sources in the sky 
should therefore exhibit an approximate large-scale dipolar distribution pointing 
towards the same direction. 

Using the AGNs in the Piccinotti et al (1982) flux-limited sample of X-ray 
sources (2-10 kcV flux limit - 3 x 10" 11 erg cm" 2 s" 1 ), Miyaji & Boldt (1990) and 
Miyaji (1994) found the dipole of these sources to point towards / = 318°, b = 
38° with a large error circle (~ 30° radius). The dipole appears to saturate at 
50—100 hr 1 Mpc and is roughly aligned with the CMB dipole. Within the framework 
of linear theory, this allows the bias parameter of the X-ray selected AGN to be 
estimated, giving a somewhat large value (bx^Q ' 6 ~ 3 — 6). Uncertainties come 
primarily from the indctcrmination of the rcdshift at which the dipole saturates. 

Plionis & Kolokotronis (1998) and Kolokotronis et al (1998) have measured the 
dipole of an X-ray flux-limited sample of galaxy clusters. This is again in rough 
alignment with the CMB dipole, but it appears to saturate at ~ 160/i _1 Mpc. 
As expected in all popular scenarios where clusters arise in extreme peaks of the 
underlying dark-matter distribution, they exhibit a large bias parameter (bx ~ 4, 
see Table 4). 

The fact that the dipoles of the two most numerous classes of extragalactic X-ray 
sources (AGNs and clusters) are roughly aligned with the CMB dipole is encour- 
aging. We note, however, that all-sky deeper samples of these objects (particularly 
X-ray selected AGN) would enormously help in defining the distance at which the 
contribution to the dipole saturates and therefore in measuring the bias parameter. 

2.2. The dipole of the X-ray background 

There are two reasons why the XRB should show a dipole signal: our motion rela- 
tive to the CMB rest frame (the so-called Compton-Getting effect) and the excess 
contribution of the sources that cause this motion in the same direction. The XRB 
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dipole is expected to be aligned with the CMB dipole, but the amplitude should be 
larger than the Compton-Getting effect, allowing for the excess emissivity. 

There are two basic problems in measuring the XRB dipole: one is the contribu- 
tion of the Galaxy and the other one is the integrated nature of the XRB whereby 
confusion noise dominates on all angular scales. Warwick, Pye & Fabian (1980) 
realized that even at photon energies > 2 keV and galactic latitudes b |> 20° a 
residual galactic contribution ~ 2 — 7% is present. Iwan et al (1982) modelled this 
galactic component in terms of a finite radius disk with thermal spectrum at T ~ 9 
keV. To emphasize how difficult is to obtain the extragalactic signal, the galactic 
contribution amounts to a few % at the galactic poles, while the effect it is being 
looked for is less than 1%. 

Attempts to look for singular enhancements of the XRB surface brightness in- 
clude those by Warwick et al (1980), the Jahoda & Mushotzky (1989) search for 
emissivity from the great attractor, the Mushotzky & Jahoda (1992) search for 
XRB negative fluctuations towards the most prominent voids and the unsuccessful 
detection of X-ray emission from superclusters by Persic et al (1990). 

By modelling out the Galaxy, Shafer (1983) and Shafer & Fabian (1983) found 
a dipole signal significant at <~ 2a level in the HEAO-1 A2 map. Most of the sub- 
sequent dipole refinements have used the same data with increasingly finer correc- 
tions for detector drifts and other unwanted effects. The latest one is by Scharf et 
al (1999), who excluded the galactic plane, the Magellanic clouds and also regions 
around the Piccinotti et al (1982) sources, which leaves less than 50% of the sky 
for the dipole analysis. Various methods are used to deal with the masked regions 
(including spherical harmonic reconstruction) and the results are shown in Table 
1. The dipole signal is very clearly detected and its intensity appears larger than 
the Compton-Getting effect. The direction of this extra large-scale structure dipole 
caused by the fluctuations in the source density is only roughly aligned with the 
direction of our motion, and its amplitude is similar to the Compton-Getting effect 
as predicted by theory (Lahav, Piran & Treyer 1997). 

In an analysis of the ROSAT all-sky data (0.9-2.4 keV), Plionis & Georgan- 
topoulos (1999) also find a dipole component. The Galaxy is modeled according to 
the Iwan et al (1982) model and they further exclude other regions associated with 
the Galaxy. The direction of the resulting dipole is in better agreement with the 
CMB dipole, but the amplitude is almost a factor of 10 larger than the Compton- 
Getting effect. 

There are various reasons for the discrepancy between these measurements. 
First, an extra residual contribution from the Galaxy is likely to contaminate more 
strongly the ROSAT data than the HEAO-1 A2 data. This would explain why the 
ROSAT dipole points closer to the galactic plane and that its amplitude is larger. A 
second reason for the discrepancy is the fact that Scharf et al (1999) have excluded 
regions around the galaxy clusters present in the Piccinotti et al (1982) sample 
(which arc known to have a very large bias parameter and represent 50% of the 
extragalactic sources in that sample) but Plionis & Georgantopoulos (1999) have 
not. In fact these last authors note that the contribution from the Virgo cluster 
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alone is of the order of 20% of the detected dipole. A good exercise which could give 
some insight on the level of the galactic contamination in the ROSAT data would 
be to exclude the clusters in the ROSAT analysis and not excise the Piccinotti ct 
al (1982) sources from the analysis of the A2 data. 

2.3. Higher order multipoles of the X-ray background 

Lahav, Piran & Treyer (1997) proposed the use of a multipole expansion of the 
angular variations of the XRB in order to measure the large-scale structure of the 
Universe. Under fairly general assumptions, the coefficients ai m of the harmonic 
expansion would be the sum of a large-scale structure term a]^^ oc l~ aA and a 
confusion noise term which is a function of the flux S cu t down to which sources have 
been excised from the maps for the multipole analysis a^ 01 '^ oc S],^ 1 , where 7 is 
the slope of the integral source counts in the energy band used (N(> S) oc 5~ 7 ). 

Treyer et al (1998) performed this analysis on the HEAO-1 A2 all sky data by 
removing regions around the Piccinotti et al (1982) sample and the galactic plane. 
They find evidence for a growth of the spherical harmonic coefficients growing at 
low values of I in a manner roughly consistent with the predictions. The significance 
of the signal is difficult to assess as the harmonic coefficients are not independent 
due to cross-talk between different orders introduced by the masking. Assuming a 
rcdshift dependent bias parameter for the X-ray sources parametrized as bx(z) = 
bx(0) + z[bx(0) — 1] (which assumes that all galaxies form at some past epoch, Fry 
1996), they estimate a rather modest bias parameter (1.0 < bx(0) < 1-6). In their 
diagrams it is also seen that the dipole (1 = 1) has an unusually large amplitude 
compared to higher harmonics. 

The way to go is indeed to have precise measurements of the XRB intensity 
on large angular scales, but with the possibility of excluding sources down to the 
faintest possible levels. Treyer et al (1998) suggest that an all-sky map with XRB 
intensities measured with a 1% precision and with sources excised down to 3 x 
10~ 13 ergcm~ 2 s^ 1 (i.e., 100 times fainter than the Piccinotti et al catalogue) would 
be ideal for a spherical harmonic analysis. 

3. CROSS-CORRELATIONS OF GALAXY CATALOGUES WITH XRB INTEN- 
SITIES 

An alternative way that has been devised to look for structure in the X-ray sky 
is to cross-correlate the unresolved XRB intensity with catalogues of galaxies. The 
amplitude of the cross-correlation function (CCF) between the X-ray intensity Ixrb 
and the galaxy surface density N g , W Xg (0) = {IxRBN g ) g / (Ixrb){N 9 ) - 1 at zero- 
lag (6 = 0) provides an approximate measurement of the fraction of the XRB arising 
either in the catalogued galaxies or in sources clustered with them within a scale of 
the beam with which X-ray observations have been obtained (Lahav et al 1993). 

Positive signals have been found for Wx g , typically of the order of 1% when the 
galaxies are optically or infrared selected, and up to > 10% when active galaxies 
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XRB data 


Galaxy catalogue 


Wxg(O) 


Rcl 


HEAO-1 A2 


IRAS 2Jy 


7 x KT 3 


M94 


GINGA NGP 


IRAS 0.7Jy 


1.4 x KT 2 


C95 


GINGA NGP 


UGC 


1.1 x KT 2 


C95 


HEAO-1 A2 


IRAS 12/z all 


< 9.6 x lO- 3 


B95 


HEAO-1 A2 


IRAS 12/i Syl 


1.1 x KT 1 


B95 


HEAO-1 A2 


IRAS 12/i Sy2 


3.1 x KT 2 


B95 



TABLE 2. Cross-correlation signal at zero-lag {Wx g ) of galaxy catalogues with 
XRB data either from the HEAO-1 A2 experiment or the GINGA North Galactic 
Pole scans. References are M94: Miyaji et al (1994); C95: Carrera et al 1995; (B95): 
Barcons et al (1995). 



are selected (see table 2). The interpretation of this signal requires to model the 
clustering of X-ray sources around the catalogued galaxies, which is indeed mod- 
ulated by the bias parameter bx- Using bx = 1, it is found that the local volume 
emissivity of optically selected galaxies amounts to ~ 10 39 Ziergs -1 Mpc~ (Lahav 
et al 1993; Miyaji et al 1994; Carrera et al 1995), most of which is contributed by 
Seyfert galaxies and QSOs (Barcons et al 1995). 

When this volume emissivity is extrapolated to higher rcdshifts, the fraction 
of the XRB intensity due to the precursors of the catalogued galaxies can be pre- 
dicted (Lahav et al 1993). Carrera et al (1995) find that 10-30% of the hard X-ray 
background might be produced by optically selected galaxies without exceeding the 
upper limits on the autocorrelation function of the XRB. This value is similar to the 
result of cross-correlation analyses of deep ROSAT X-ray images with deep optical 
images in the same fields (Almaini et al 1997). 

A constraint on the bias parameter of X-ray sources from the CCF results can be 
derived by taking into account that a fraction / w 2/3 of the CCF signal arises from 
sources clustered with the catalogued galaxies. As that contribution scales linearly 
with bx, the local volume emissivity scales oc 1+ ^ bx ■ Since the AGN-only local 
emissivity is also ~ 10 39 ergcm~ 2 s _1 , we can safely derive that bx < 2 as otherwise 
the total volume emissivity will be significantly less than the AGN emissivity. 

4. CLUSTERING OF X-RAY SELECTED SOURCES 

In the recent years large, complete samples of X-ray selected AGN have been built. 
This has allowed the direct measurement of the 3D spatial correlation function £ (r) 
for these objects and its comparison with the spatial correlation function of galaxies 
selected at other wavebands. 

Carrera et al (1998) used two complete samples of X-ray selected AGN in pencil 
beam survey regions, spanning a wide redshift range (0 < z < 2) to search for 
clustering signals and deriving its amplitude and redshift evolution. Clustering is 
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found to be 99% significant at z < 1. When the spatial correlation function is fitted 
to a standard power-law form £(r) = (1 + 2) p (^) -1 ' 8 (for comoving r), it is seen 
that comoving or slower clustering evolution is excluded, and that even for stable 
or linear growth the values of tq permitted by the data are of the same order as the 
ones derived from clustering of IRAS galaxies. Carrera et al (1998) conclude that 
X-ray selected AGN are not significantly biased 0.7 < bx < 2. 

Akylas, Plionis & Georgantopoulos (1999) have used the ROSAT All Sky Survey 
sources to derive a local angular correlation function from which they estimate a 
somewhat higher correlation length (ro ~ 7 — 9Mpc), consistent with optically 
selected QSO clustering (La Franca et al 1998) and comoving clustering evolution. 
The obvious weakness of this method is that it is not based on 3D but 2D data. 

5. FLUCTUATIONS AND ANISOTROPIES IN THE XRB 

The method of auto-correlating the XRB intensity at various separations has been 
extensively used in an effort to detect small scale structure in the XRB attributable 
to source clustering (Barcons & Fabian 1989, De Zotti et al 1990, Jahoda & Mushotzky 
1991, Carrera et al 1991, Carrera & Barcons 1992, Carrera et al 1993, Chen et al 
1994). These works produced a set of upper limits for the auto-correlation function of 
the XRB Wxx{&) — {IxrbIxrb)sI '(Ixrb) 2 — 1 on different angular scales (except 
for the Jahoda & Mushotzky 1991 work, which claimed a detection at separations 
~ 10°) of the order of 10~ 3 — 10~ 4 which constrained the clustering properties of 
the underlying source population (see, e.g., Fabian & Barcons 1992). 

Under the assumption of comoving clustering evolution, the sources of the XRB 
cannot be more strongly clustered than optically selected galaxies (see Carrera & 
Barcons 1992), in which case bx ~ 1- However, as explained above, Carrera et al 
(1998) found marginal evidence for faster clustering evolution in samples of X-ray 
selected AGN. This means that b x could be higher without violating the upper 
limits on the autocorrelation function, as the sources that produce the bulk of the 
XRB at high redshift could be very weakly clustered. 

A further method employing the XRB angular variations has been to search 
for fluctuations in the XRB intensity distribution in excess of the ones expected 
from confusion noise produced by unresolved sources. These excess fluctuations 
should then be attributed to source clustering if all remaining noises (counting 
noise, systematics, etc.) could be removed. Studies of this kind have invariably 
lead to upper limits summarized in Table 3. What actually limits the sensitivity 

— 1/2 

of this method is the statistics: it scales as N ob J where N bs is the number of 
independent measurements of the XRB intensity that have been used to derive the 
excess fluctuations. 

Excess fluctuations are related to the power spectrum of the density field of 
the Universe, weighted with the X-ray volume emissivity as a function of redshift 
(Barcons, Fabian & Carrera 1998). The method is potentially very powerful as it 
reflects the clustering properties of the sources that produce the bulk of the XRB 
at redshifts z > 1. 
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XRB data 


Beam 


(M) 

\ I I races.! 


Rcf 


HEAO-1 A2 


5° x 5° 


< 0.02 


S83 


Ginga 


2° x 1° 


< 0.04 


B97 


ROSAT 


7T x (2.5') 2 


< 0.07 


CFB97 


ROSAT 


^(10' 2 - 5' 2 ) 


< 0.12 


CFB97 



TABLE 3. Upper limits to excess fluctuations. References are S83: Shafer (1983); 
B97: Butcher et al (1997); CFB97: Carrera, Fabian & Barcons (1997) 



Measurement 


Reference 


Scale (Mpc) 


b (fi = 1) 


X-ray cluster dipole 


PK98 


10-100 


4 


X-ray AGN dipole 


MB90 


1000 


3-6 


XRB-galaxy CCF 




10-100 


< 2 


XRB dipole vs bulk motions 


S99 


1000 


2-7 


XRB multipoles vs bulk motions 


T98 


100-1000 


1-2 


Clustering of distant AGN 


C98 


10-100 


1-2 


Clustering of nearby AGN 


APG99 


100 


2-3 



TABLE 4. Bias parameters as inferred from various measurements, references are: 
PK98 Plionis & Kolokotronis (1998); MB90: Miyaji & Boldt (1990); S99: Scharf 
et al (1999); T98: Treyer et al (1998); C98: Carrera et al (1998); APG99: Akylas, 
Plionis & Georgantopoulos (1999). 

6. THE BIAS PARAMETER OF X-RAY SOURCES 

Over the past sections we have discussed various approaches to detect and measure 
the clustering properties of X-ray sources. Tabic 4 summarizes the inferred bias 
parameter bx from these studies. Measurements are carried out with a variety of 
methods, correspond to different objects, are sensitive to different redshifts and 
also to different scales. Besides that, all dynamical estimates actually measure the 
combination bx^ ' 6 - 

Measurements of the correlation function are also affected by the cosmological 
parameters in the computation of the distances at significant redshifts, beyond the 
obvious linear dependence on Hq. If we live in an accelerating Universe, the Carrera 
et al (1998) correlation length would have to be scaled up by 30-50%, resulting 
in a subsequent increase of almost a factor of 2 in the bias parameter. Given the 
uncertainties in the values of qo and A (even for a flat Universe), the Carrera et al 
(1998) and Akylas et al (1999) results cannot be considered inconsistent. 

As expected, clusters are a largely biased population (bx ~ 4) compared to AGN 
(bx ~ 1 — 2). The multipoles of the XRB are expected to be dominated by AGN, as 
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these objects are the main sources of the XRB. The bias parameter derived from the 
XRB multipolcs is consistently in agreement with the bias parameter derived from 
AGN clustering (bx ~ 1 — 2). The exception to this is the XRB dipolc which implies 
a larger value of bx- This could be partly due to a larger cluster contribution, as 
the lowest order multipoles are most sensitive to nearest (and brightest) sources, 
where the cluster contribution to the source counts (~ 10% on average in the deep 
extragalactic surveys) is ~ 50% for the Piccinotti et al (1982) sample. 

7. FUTURE PROSPECTS 

X-ray astronomy is now in a position to address cosmological studies. X-ray selected 
AGN which produce most of the X-rays in the Universe, appear to trace mass 
with a moderate bias parameter bx ~ 1 — 2, but that has to be better defined as 
a function of scale and redshift. Chandra and XMM will carry out several deep 
'pencil beam' surveys which, after subsequent identification of the serendipitous 
sources discovered, will define the redshift evolution of the AGN X-ray luminosity 
function at photon energies > 2keV and therefore the X-ray volume emissivity as 
a function of redshift. However, these surveys will not map sufficiently large areas 
of the sky which are necessary to trace the large-scale structure of the Universe at 
the redshifts where the XRB was produced. 

The obvious way to go would be to survey very large areas of the sky (the whole 
sky even better) for X-ray sources, in order to have a most complete picture. Unless 
hard X-rays are produced at significantly lower redshifts than soft X-rays (which 
is doubtful in view of the ASCA and BeppoSAX surveys), to reach z ~ 1 where 
a significant fraction of the X-ray emissivity in the Universe resides, these surveys 
will have to go at least down to ~ 10~ 14 crgcm~ 2 s _1 . 

There is an alternative which is to perform high sensitivity observations of the 
XRB with a beam corresponding to the linear scale to be probed (Barcons, Fabian & 
Carrera 1998). As the peak of the power spectrum of the density field of the Universe 
occurs at comoving wavenumbers ~ 0.01 — 0.1 hMpc^ 1 , for a standard geometry a 
1° resolution is well matched to this at z ~ 1 — 3. All-sky measurements of the XRB 
intensity on that angular scale with a precision of a few % could then be used to 
detect the excess fluctuations due to source clustering which are expected to be just 
below 1% in amplitude. Controlling all other possible sources of excess fluctuations 
well below that level requires a stable large-area detector (to reduce photon counting 
noise) and probably an X-ray monitor which images simultaneously the brightest 
sources in the field. 
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